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In this work, we have reviewed the Oslo method, which enables the simultaneous extraction of level density 
and y-ray transmission coefficient from a set of particle-y coincidence data. Possible errors and uncertainties 
have been investigated. Typical data sets from various mass regions as well as simulated data have been tested 
against the assumptions behind the data analysis. 
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I. INTRODUCTION 

Nuclear level densities and 7-ray strength functions are two 
indispensable quantities for many nuclear structure studies 
and applications. The approach developed by the nuclear 
physics group at the Oslo Cyclotron Laboratory (OCL), was 
first described in 1983 [1 1, and has during a peiiod of nearly 
30 years been refined and extended to the sophisticated level 
presently referred to as the Oslo method. The method is de- 
signed to extract the nuclear level density and the 7-ray trans- 
mission coefficient (or 7-ray strength function) up to the neu- 
tron (proton) threshold from particle-7 coincidence data. Typ- 
ical reactions that have been utilized are transfer reactions 
such as (^He, ay) and (p,ty), and inelastic scattering reac- 
tions, e.g., (^He, ^He'7) and (p,p'7). 

The measurements have been very successful and have shed 
new light on important issues in nuclear structure, such as 
the Ml scissors mode lIlKSl, and the sequential breaking of 
Cooper pairs ||6] |7]. For the Sn isotopes, a resonance-like 
structure that may be due to the so-called E\ pygmy reso- 
nance has been observed JH ID. Also, a new, unpredicted 
low-energy increase in the y-ray strength function for Ey <3 
MeV of medium-mass nuclei has been discovered by use of 
this method |10-15|. This enhancement might have a non- 
negligible impact on stellar reaction rates relevant for the nu- 
cleosynthesis Iil6j . At present, this structure is poorly under- 
stood. 

In this work, we have investigated the possible systematic 
errors that can occur in the Oslo method due to experimen- 
tal limitations and the assumptions made in the various steps 
of the method. We have studied experimental data with high 
statistics (and thus low statistical errors) so that systematic er- 
rors can be revealed. Also, we have used simulated data to en- 
able better control on the input parameters (level density and 
7-ray strength function). In Sec. |Il]we give a shoit overview 
of the experiments and the various steps in the method. In 



Sec. Ill we present the possible systematic eiTors for each 
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main step of the method. Finally, a summary and conclud- 
ing remarks are given in Sec.|IV| 



II. EXPERIMENTAL PROCEDURE AND THE OSLO 
METHOD 

The Oslo method is in fact a set of methods and analysis 
techniques, which together make it possible to measure level 
density and 7-ray strength from particle-7coincidence data. In 
this section, these techniques and methods wiU be described. 



A. Experimental details 

The experiments were performed at the OCL using a light- 
ion beam delivered by the MC-35 Scanditronix cyclotron. 
Typically, ^^He beams with energy 30 — 45 MeV has been used. 
Recently, also proton beams with energy 15 — 32 MeV have 
been applied. Self-supporting targets enriched to w 95% in the 
isotope of interest, and with a thickness of sa 2 mg/cm^ were 
placed in the center of the multi-detector array CACTUS 1 17 |. 
CACTUS consists of 28 colUmated Nal(TI) 7-ray detectors 
with a total efficiency of 15.2(1)% for By = 1332 keV. Usu- 
ally also a 60% Ge detector has been applied to monitor the 
populated spin range (^ — %K) and possible target contam- 
inations. The experiments were typically run for a period of 
1 —2 weeks with beam currents of ^ 1 nA. 

Inside the CACTUS array, eight collimated Si particle de- 
tectors were used for detecting the charged ejectiles from the 
nuclear reactions. The particle detectors were placed at 45° 
relative to the beam line in forward direction. The detectors 
were of — E type with a thin (^ 140 pLm) front detec- 
tor and a thick (^ 1500 /im) end detector where the charged 
ejectiles stop. The particle telescopes enable a good separa- 
tion between the various charged-ejectile species. The energy 
resolution of the particle spectra ranges from 150 — 300 keV 
depending on the beam species, the mass of the target nucleus 
and the size of the collimators. Both singles and coincidence 
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FIG. 1: Singles a-particle spectrum (upper panel), and a-y coinci- 
dence spectrum (lower panel) from the ' '^Sn(-'He, o;)"^Sn reaction. 
The data are taken from the experiment presented in Ref. |7|. 



FIG. 2: Singles He spectrum (upper panel), and He-y coincidence 
spectrum (lower panel) from the '^'*Dy(-'He,-'He')'^'*Dy reaction. 
The data are taken from the experiment presented in Ref. IS]. 



events were measured for the ejectiles. In Figs. [T| [2] and [3] 
the singles and particle-7 coincidence spectra are shown for 
the reactions "^Sn(3He,a)^i'^Sn, ^^'^DyClle,^He'y^'^Dy and 
^^Ti(p, p')'^^Ti, respectively. 

The ejectile energy can easily be transformed to the excita- 
tion energy of the residual nucleus using the reaction Q-values 
and kinematics. Thus, an excitation energy vs. 7-ray energy 
matrix can be built from the coincidence events. An exam- 
ple of such a matrix is shown in Fig. |4] (left panel), where 
the 7-ray spectra have been corrected for the known response 
functions of the CACTUS array 1 1 8 1 . The correction (or un- 
folding) method is described in detail in Ref. ifTSll . One of 
the main advantages with this method is that the fluctuations 
of the original spectra are preserved without introducing addi- 
tional, spurious fluctuations. 

B. Extracting primary y rays 

As 7 decay from highly excited states often involves a cas- 
cade of transitions, it is necessary to isolate the 7 rays that are 
emitted in the first decay step of all the possible decay routes, 
since information on the level density and the 7-ray strength 
function can be extracted from the distribution of these pri- 
mary 7 rays (also called first-generation 7 rays). Therefore, a 
method has been developed in order to separate the primary 
7-ray spectra from the 7 rays that origin from the later steps 
in the decay cascades at each excitation energy. This method 



will hereafter be referred to as \he first-generation method and 
is described in detail in Ref. ||T9| . This method is very impor- 
tant, since the correctness of the further analysis is completely 
dependent on correctly determined primary 7-ray spectra. The 
first-generation method shares several features with another 
subtraction technique developed by Bartholomew et al. |20i . 
The main features of the method will be outlined in the fol- 
lowing. 

From the E vs. Ey matrix, where 7-ray spectra for each 
excitation energy are contained, the primary 7-ray spectra are 
extracted through an iterative subtraction technique. The un- 
folded spectra fi are made of all generations of 7 rays from all 
possible cascades decaying from the excited levels within the 
excitation-energy bin /. Now, we utilize the fact that the spec- 
tra fj^i for all the underlying energy bins j contain the same 
7- transitions as /, except the first 7 rays emitted ^ since they 
will bring the nucleus from the states in energy bin / to under- 
lying states in the bins j. Thus, the primary 7-ray spectrum hi 
for each bin / can be found by 

hi=fi-gi, (1) 



This is only true if tlie y-decay pattern is the same regardless of whether 
the states in the bin were populated from the direct reaction, or from decay 
from above-lying bins. See SeclIIIBI 
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FIG. 3: Singles proton spectrum (upper panel) and p-y coincidence 
spectrum (lower panel) from the '*''Ti(p, p')'*^Ti reaction. The data 
are taken from the experiment presented in Ref. L15I . 



where is a weighted sum of all spectra 
gi = nnwnfi + naWiifi + ■■.+ riijWijfj = ^ ntjWijfj . (2) 

Here, the unknown coefficients Wij (with the normalization 
Y^j^ij = 1) represent the probability of the decay from states 
in bin / to states in bin j. In other words, the Wij values make 
up the weighting function for bin /, and contain the distribu- 
tion of branching ratios as a function of 7-ray energy. There- 
fore, the weighting function w, corresponds directly to the pri- 
mary 7-ray spectrum /i, for bin /. 

The coefficients rtij are correcting factors for the different 
cross sections of populating levels in bin / and the underlying 
levels in bin j. They are determined so that the total area of 
each spectrum /, multiplied by n,y corresponds to the same 
number of cascades. This can be done in two ways |fT9l : 

• Singles normalization. The singles-particle cross sec- 
tion is proportional to the number of events populating 
levels in a specific bin, and thus to the number of de- 
cay cascades from this bin. We denote the number of 
counts measured for bin / and j in the singles spectrum 
Si and Sj, respectively. The normalization factor n,j that 



should be applied to the spectrum fj is then given by 



Si' 



(3) 



Assume an A^-fold population of an excited level E. The 
decay from this level will result in 7-ray cascades, 
where the ith cascade contains M, 7 rays. The average 
7-ray energy (Ey) is equal to the total energy canied by 
the 7 rays divided by the total number of 7 rays: 

E EE 

Wr 



-N- 



Ll.Mi 



(4) 



Then, the average 7-ray multiplicity is simply given by 
E 



{M) = 



(5) 



"7/ 



The average 7-ray multiplicity (M,) can thus easily be 
calculated for each excitation-energy bin /. Let the area 
(or total number of counts) of the 7-ray spectrum /, be 
denoted by A{fi). Then the singles particle cross sec- 
tion Si is proportional to the ratio A{ fi)/{Mi), and the 
normalization coefficient riij that should be applied to 
bin / when subtracting bin j is 

A{fi)/{M.) _ {Mj)A{fi) 



(6) 



A{fj)/{Mj) {M.)A{fj)- 
The two normalization methods give normally the sa me re 
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suits within the experimental error bars-, see also Sec 
The resulting primary 7-ray matrix of ^"V is shown in Fig. [4] 
(right panel), using the singles normalization method. 

In cases where the multiplicity is well determined, an area 
consistency check can be applied to Eq. ([T]l. Assume that a 
small coiTection is introduced by substituting gi by 8gi, where 
8 is close to unity. The area of the first-generation 7 spectrum 
is then 



A{h,)=A{fi)-5A{gi), 



(7) 



and this coiTesponds to a 7-ray multiplicity of one unit. Since 
the number of primary 7 rays in the spectrum fi equals 
A{fi)/{Mi),A{hi) is also given by 



Aihi)=A{fi)/{M.). 
Combining Eqs. (|7]l and ^ yields 

5 = (1-1/(M,))^. 



(8) 



(9) 



The 5 parameter can be varied to get the best agreement of the 
areas A(/i,), A{fi) and A{gi) within the following restriction: 
5 = 1.00±0.15; thatis, the correction should not exceed 15%. 
If a larger coiTection is necessary, then improved weighting 
functions should be determined instead. 

As mentioned before, the weighting coefficients Wij cor- 
respond directly to the first-generation spectrum /i,, and this 
close relationship makes it possible to determine Wij (and thus 
hi) through a fast converging iteration procedure lfT9ll : 



• Multiplicity normalization. The average 7-ray multi- 
plicity (M) can be obtained in the following way II2TI : 



^ In case of the presence of isomeric states, the multiplicity method must be 
used to get the correct normahzation. 
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FIG. 4: (Color online). Excitation-energy vs. y-ray energy matrix (left), and primary y-ray matrix (right) for V. The data are taken from the 
experiment presented in Ref. 1131 . 



1. Apply a trial function for w,y. 

2. Deduce hi. 

3. Transform hi to by giving /i, the same energy cali- 
bration as Wij, and normalize the area of hi to unity. 

4. If w,y(new) sa w,j(old), convergence is reached, and the 
procedure is finished. Otherwise repeat from step 2. 

The first trial function could be the unfolded spectrum //, or a 
theoretical estimate based on a model for the level density and 
7-ray transmission coefficient, or a constant function; it turns 
out that the resulting first-generation spectra are not sensitive 
to the starting trial function. Also, previous tests of the conver- 
gence properties of the procedure have shown that excellent 
agreement is achieved between the exact solution (from sim- 
ulated spectra) and the trial function Wij already after three it- 
erations |[T9l . Usually, about 10 — 20 iterations are performed 
on experimental spectra. 



C. Determining level density and y-ray strength 

Once the primary y-ray spectra are obtained for each exci- 
tation energy, the first-generation matrix PiE^Ey) is used for 
the determination of level density and y-ray strength. For sta- 
tistical y-decay, the decay probability (given by P{E,Ey)) of 
a y-ray with energy Ey decaying from a specific excitation en- 
ergy E is proportional to the level density p (Ef) at the final 
excitation energy Ef = E— Ey, and the y-ray transmission co- 
efficient ^{Ey): 



P{E,Ey)ocp{Ef)^{Ey 



(10) 



The above relation holds for decay from compound states, 
which means that the relative probability for decay into any 



specific set of final states is independent on how the com- 
pound nucleus was formed. Thus, the nuclear reaction can 
be described as a two-stage process, where a compound state 
is first formed before it decays in a manner that is independent 
of the mode of formation [.22. ,23 J . This is believed to be ful- 
filled at high excitation energy, even though direct reactions 
are used, as already discussed previously. Equation 
also be compared with Fermi's golden rule: 



10 1 can 
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(f|//,„t|i)|V(£f), 



(11) 



where A is the decay rate of the initial state |i) to the final 



state |f), and i^int is the transition operator. In Eq. ( 10 1, an en 



semble of initial and final states within each excitation-energy 
bin is considered, and thus we obtain here the average decay 
properties of a set of initial states to a set of final states. Note, 
however, that in contrast to Fermi's golden rule where the ma- 
trix element is stiictly dependent on the initial and final state, 
the transmission coefficient ^ is only dependent on the y-ray 
energy and neither the initial nor the final excitation energy. 
This is in accordance with the Brink hypothesis |24|, which 
states that the collective giant dipole mode built on excited 
states has the same properties as if built on the ground state. 
In its generalized form, this hypothesis includes all types of 
collective decay modes. Assuming that this hypothesis holds, 
the first-generation matrix P{E,Ey) is separable into two func- 
tions p and 3^ as given in Eq. ([TO]). 

To extract the level density and the y-ray transmission co- 
efficient, an iterative procedure ll25l is applied to the first- 
generation matrix P{E,Ey). The basic idea of this method 
is to minimize 
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f P,k{E,Ey)-P{E,Ey) 
V AP{E,Ey) 



(12) 



where A^fiee is the number of degrees of freedom, and 



5 



M'{E,Ey) is the uncertainty in the experimental first- 
generation 7-ray matrix. The fitted first-generation 7-ray ma- 
trix Pi-ti{E,Ey) can theoretically be approximated by 



p{E-Ey).^{Ey 



,p{E-Ey)^{Ey) 



(13) 



the maximum value Zs™^" = £ at this excitation-energy bin 
is unity: 



the p function is adjusted to fit the number of known discrete 
levels at low excitation energy and neutron (or proton) reso- 
nance data at high excitation energy. This normalization is 
shown for '^^Dy in Fig. [s] The data point at high excitation 
energy (open square in Fig. [S} is calculated in the following 
way according to Ref. ||25| . The starting point is Eqs. (4) and 
(5) of Ref. 123: 



The experimental matrix of first-generation 7 rays is normal- 
ized II25I such that for every excitation-energy bin E, the sum 
over all 7 energies Ey from some minimum value EIJ"" to Pi^T^) 



0rexp(2v^) (27 + l)exp [-(/+ l/2)V2cj2 



and 
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(17) 



£ P{E,Ey)^l. 



(14) 
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P{U) = 



fn exp (2 
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(18) 



The experimental matrix P{E.Ey) and the fitted matrix 
P(n{E,Ey) of ^^Ti are displayed in Fig.js] The energy limits 
set in the first-generation matrix for extraction are also shown. 
These limits (E™'", E™^", and Z?™") are chosen to ensure that 
the data utilized are from the statistical excitation-energy re- 
gion and that no 7-lines stemming from, e.g., yrast transitions, 
are used in the further analysis. Note that the 7-ray energy bins 
are now re-binned to the same size (120 keV in this case) as 
the excitation-energy bins. 

Each point of the p and ^ functions is assumed to be an in- 



dependent variable, so that the reduced % of Eq. ( 12 1 is min 



imized for every argument E — Ey and E. The quality of the 
procedure when applied to ^^Ti and '^^Dy is shown in Figs.|6] 
and|7] where the experimental first-generation spectra for var- 
ious initial excitation energies are compared to the least-;)^^ 
solution. In general, the agreement between the experimental 
data and the fit is very good. Note, however, that differences 
of several standard deviations do occur In the case of ^^Ti, 
this is particularly pronounced for the peaks at Ey sa 4.7 MeV 
for E = 5.6 MeV and Ey w 5.5 MeV for E = 6.4 MeV. As 
these peaks correspond to the decay to the first excited state, 
one might expect large Porter- Thomas fluctuations ll26l in the 
strength of these transitions. 

The global fitting to the data points only gives the functional 
form of p and In fact, it has been shown f25\ that if one so- 
lution for the multiplicative functions p and is known, one 
may construct an infinite number of other functions, which 
give identical fits to the P{E,Ey) matrix by 



p{E-Ey) = At\.p[a{E-Ey)]p{E-E. 
.7{Ey) = Bexp(a£'y)^(£'y). 



(15) 
(16) 



Therefore the transformation parameters a, A and B, which 
correspond to the physical solution, must be found from ex- 
ternal data. 



D. Normalization 

In order to determine the correction a to the slope of the 
level density and the 7-ray transmission coefficient, and to de- 



where p (t/,7) is the level density for a given spin J , and p (t/) 
is the level density for all spins. The intrinsic excitation en- 
ergy JJ , the level-density parameter a, and the spin cutoff 
parameter a are normally taken from Ref. 1281 in previous 
works, or from Ref. l29l in recent works. 

Now, let 7, be the spin of the target nucleus in a neutron res- 
onance experiment. The average neutron resonance spacing 
D^'=o for s-wave neutrons can be written as 



1 



1 



[p(5„,y = /, + l/2)+p(5„,7 = /,-l/2)], (19) 



because all levels with spin 7 = ± 1/2 are accessible in neu- 
tron resonance experiments, and because it is assumed that 
both parities contribute equally to the level density at the neu- 
tron separation energy 5„. Combining Eqs. (17i-(19i with 
IJ = S„, one finds the total level density at the neutron separa- 
tion energy to be 



piS„ 



2a'- 



1 



Do {I, + l)exp[-(/, + l)2/2ff2] +/,exp [-1? /la^] 

(20) 

Note also that the resonance spacing between p-waves, Di, 
could also be used for calculating p{S„), see Ref. |30|. 

Since our experimental data only reach up to excitation en- 
ergies around S„ — /s™", an interpolation has been made be- 
tween the Oslo data and p(5„) using the back-shifted Fermi 
gas model of Refs. l28l |29l , as shown in Fig. [8] It should 
be noted that in most cases the gap between the data and 
p{S„) is small, so that the normalization is not very sensi- 
tive to the interpolation; a pure exponential function of the 
type p{E) — Coexp(Ci£'), where Co and C\ are fitting param- 
eters, gives an interpolation of equally good agreement (see 
Ref. 0). 

The slope of the 7-ray transmission coefficient S/'iEy) has 
akeady been determined through the normalization of the 
level density, as explained above. The remaining constant B 
in Eq. ( 16 1 gives the absolute normalization of 5^, and it is 



determined using information from neutron resonance decay 



termine the absolute value A of the level density in Eq. ( 15 1, 



on the average total radiative width 
Ref. ED. 



Ty) at 5„ according to 
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FIG. 5: (Color online). Experimental first-generation matrix P{E,Ey) (left) and the calculated P±{E,Ey) (right) of ^^Ti from the iteration 
procedure of A. Schiller et al. 1251 . The dashed lines show the limits set in the experimental first-generation matrix for the fitting procedure. 
The data are taken from the experiment presented in Ref. 1151 . 
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FIG. 6: (Color online). Experimental first-generation y spectra (data points with error bars) at six different initial excitation energies (indicated 
in the figure) compared to the );} fit (solid lines) for "^^Ti. The fit is performed simultaneously on the entire first-generation matrix of which 
the six displayed spectra are a fraction. The data are taken from the experiment presented in Ref. II5I . 



The starting point is Eq. (3.1) of Ref. 

{Yy{E,j,K))^ } ,^) iL r 

X p{E -Ey,J{,7tf), 



(21) 



where {ry{E,J, n)) is the average total radiative width of lev- 
els with energy E, spin J and parity n. The summation and 



integration are going over all final levels with spin Jf and par- 
ity Tii that are accessible through y transitions with energy 
Ey, electromagnetic character X and multipolarity L. Assum- 
ing that the main contribution to the experimental ,^ is from 
dipole radiation (L = 1), we get 



B.^{Ey) - BY^:7xL{Ey) ~ B [^Ei {Ey) + ^mi [Ey)] , (22) 
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FIG. 7: (Color online). Same as Fig.Hfor ' Dy. The data are taken from the experiment presented in Ref. \5\. 
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FIG. 8: Normalization procedure of the level density (data points) of 
I64j)y jjjg points between the arrows are normalized to known 
levels at low excitation energy (solid line) and to the level density at 
the neutron separation energy (open square) using an interpolation 
with the Fermi-gas level density (dashed line). The data are taken 
from the experiment presented in Ref. |5]. 



from which the total, experimental 7-ray strength function can 



easily be calculated: 

f{Er) 



(23) 



from the relation between 7-ray strength function and 7-ray 
transmission coefficient lf33]l: 



(24) 



^XL{Ey)^2nE^^''+'\fxL{Er). 



Further, we also assume that there are equally many acces- 
sible levels with positive and negative parity for any excitation 
energy and spin, so that the level density is given by 

p{E-Ey,Ju±%f) = ^-p{E-Ey,Ji). (25) 

Now, by combining Eqs. (j2T]l, ( [22| and ( [25] l, the average to- 
tal radiative width of neutron s-wave capture resonances with 
spins /, ± 1/2 expressed in terms of the experimental ."^ is 
obtained: 



(r^(5„,/,±l/2,7r,)) 



B 



dEy,^{Ey 



xp{S„~Ey) ^ giS„~Ey,I,±l/2+J), 

J=-\ 

(26) 

where /, and tt, are the spin and parity of the target nucleus 
in the (n, 7) reaction, and p (5„ — Ey) is the experimental level 
density. Note that the factor l/p(5„,/, ± 1/2, ;r,) equals the 
neutron resonance spacing Dq. The spin distribution of the 
level density is assumed to be given by 1271 : 



g{E,J) : 



27+1 

2C72 



exp[-(y+l/2)V2(j2] . (27) 
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The spin distribution is normalized so that ^^^(ii, 7) 1. The 
experimental value of (Fy) at S,, is then the weighted sum of 



the level widths of states with 7, ± 1/2 according to Eq. (26 1. 
From this expression the normalization constant B can be de- 
termined as described in Ref. [31 J . However, some considera- 
tions must be done before normahzing according to Eq. (|26]). 

Methodical difficulties in the primary 7-ray extraction pre- 
vent determination of the function ."^{Ey) for Ey < /i™'" as 
discussed previously. In addition, the data at the highest 7- 
energies suffer from poor statistics. Therefore, 3^ is extrap- 
olated with an exponential function, as demonstrated for ^' V 
in Fig. |9] The contribution of the extrapolation to the total 
radiative width given by Eq. ( 26 1 does not normally exceed 
15%, thus the errors due to a possibly poor extrapolation are 
of minor importance II31I . 




FIG. 9: Extrapolation of the y-ray transmission coefficient of ^'V. 
The data points between the arrows in the low and high y-energy 
regions are utilized to fit the exponential functions to the data. The 
data are taken from the experiment presented in Ref. I13J . 



III. UNCERTAINTIES AND POSSIBLE SYSTEMATIC 
ERRORS 

In the following, we will go through the Oslo method step 
by step with a close look at the uncertainties and possible sys- 
tematic eiTors connected to it. 



A. Unfolding of y-ray spectra 

The unfolding method is described in great detail in 
Ref. ifTsl . The method is based on a successive subtraction- 
iteration technique in combination with a special treatment of 
the Compton background. The stability of the method has 
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FIG. 10: (Color online). Efficiencies relative to the 1.33-MeV effi- 
ciency used in the unfolding method. 



been extensively tested in previous works and has proven to 
be very robust and reliable (see, e.g., Refs. | fT3]fT4l[30l ). To 
some extent also the impact of slightly erroneous response 
functions has been investigated in Ref. |23|. As this is the 
part of the method that has the largest potential of influencing 
the final results, it is further addressed here. 

In the unfolding method, the 7-ray spectra are corrected for 
the total absorption efficiency of the Nal crystals for a given 7 
energy. The applied efficiencies (normalized to the efficiency 
at 1.33 MeV) are shown in Fig. 



10 



We have tested the effect of reducing these efficiencies 
by up to « 20% for 7 energies above 1 MeV, using simu- 
lated particle-7 coincidences generated with the DICEBOX 
code 1 34). In the DICEBOX algorithm, a complete decay 
scheme of an artificial nucleus is generated. In this case we 
have considered an artificial nucleus resembling ^^Fe. Below 
an excitation energy of about 2.2 MeV, all information from 
the known decay scheme is used; above this energy the lev- 
els and decay properties are generated from a chosen model 
of the level density and 7-ray strength function. The code al- 
lows to take into account Porter- Thomas fluctuations of the 
individual transition intensities as well as assumed fluctua- 
tions in the actual density of levels. Each particular set of 
the level scheme and the decay intensities is called a nuclear 
realization. For more details see Ref. |34|. A restriction on 
the spin distribution of the initial excitation-energy levels of 
1/2<J< 13/2 was applied, and the chosen bin size was 120 
keV. In these simulations, each level in a bin was populated 
with the same probability independently of its spin and parity. 
This means that Porter- Thomas fluctuations were not consid- 
ered in the population of levels via the direct population, but 
only in the 7 decay. 

The results on the extracted level density and 7-ray strength 
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FIG. 12: Calculated y-ray spectra from the ^^^Dy{^lie,aY) reaction 
at 45 MeV f35 1. The solid line shows the pre-equilibrium component 
of the total y-ray spectrum (y before a), the dashed line represents 
the equilibrium part (a before y), and the filled squares give the total 
spectrum. 
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FIG. 11: (Color online). Test of sensitivity on the total y efficiency on 
the level density (top) and y-ray strength function (bottom) extracted 
from simulated data. 



function are shown in Fig. 11 It is seen that the level den 



sity is not very sensitive to the total efficiency, but the slope 
of the y-ray strength function is increased when using a too 
low efficiency for the high-energy y rays, thus leading to a 
too large coiTection in the unfolding procedure. However, it is 
very gratifying that the overall shape is indeed preserved, and 
the deviation in absolute value of the two y-ray strength func- 
tions does not exceed 20%, corresponding to the maximum 
change in the absolute efficiency. 



B. The flrst-generation method 

The first-generation method, which is applied to extract the 
distribution of primary y rays from each excitation energy, is 
a sequential subtraction technique where an iterative proce- 
dure is applied to determine the weighting coefficients 
which coiTespond to the primary y-ray spectrum as described 
in SecHTB] 

The main assumption of the first-generation method is that 
the y decay from any excitation-energy bin is independent on 
how the nucleus was excited to this bin. In other words, the 
decay routes are the same whether they were initiated directly 
by the nuclear reaction or by y decay from higher-lying states, 
giving rise to the same shape of the y spectra. This assumption 
is automatically fulfilled when states have the same cross sec- 
tion to be populated by the two processes, since y branching 
ratios are properties of the levels themselves. 

In the region of high level density, the nucleus seems to 
attain a compound-like system before emitting y rays even 
though direct reactions are utilized. This is due to two fac- 
tors. First, a significant configuration mixing of the levels will 
appear when the level spacing is comparable to the residual in- 
teraction. Second, the reaction time, and thus the time it takes 
to create a complete compound state, is « 10^'**s, while the 
typical life time of states in the quasi-continuum is w 10^ '^s. 
Therefore, it is reasonable to assume that the nucleus has ther- 
malized prior to y decay. This is supported by recent calcula- 
tions ll35l based on the Iwamoto-Harada-Bispinghoff model, 
showing that for the ^^''Dy(^He,ay) reaction with a 45-MeV 
'^He beam, the pre-equilibrium ("direct") component of the 
y-ray spectra is very small for y energies below « 10 MeV 
(see Fig. 12 1. The same is true for the '^^Ti(p,p'Y) reaction 



with proton beam energy Ep = 15 MeV, see Fig. 13 Note that 
the "direct" component is calculated using the pre-equilibrium 
(i.e. statistical) formalism. 
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FIG. 13: Same as Fig. 12 but for the inelastic scattering Ti(p,p'y) 



at£;p = 15MeV. 



Experimentally, the independence of the reaction mecha- 
nism has been tested by creating the same compound nu- 
cleus with the two different reactions (^He,a) and (^He,^He'). 
This has been done for, e.g., '^''•^^Mo lIH, i^i.ieipy g]^ ^j^^ 
171,172y^ |3l. One observes an excellent agreement with the 
level density and y-ray strength function resulting from the 
two reactions within the experimental error bars. However, at 
very low excitation energies, there is a significant difference in 
the obtained level density: the inelastic scattering gives con- 
sistently a higher level density close to the ground state than 
does the pick-up reaction. This could be a sign that the in- 
elastic scattering populates states with wave functions having 
a large overlap with the ground state and the low-lying ex- 
cited states. Thus, the decay to the ground state and low-lying 
states will be very fast, and cannot be characterized as com- 
pound decay. 

The region at low excitation energies can be tricky also 
in other aspects. Vertical ridges and/or valleys can occur in 
the primary 7-ray matrix as a consequence of differences in 
feeding of these discrete states, giving significantly different 
shapes of the 7 spectra at low E compared to higher excitation 
energies. The direct reaction cross section depends strongly 
on the intrinsic wave functions of low lying states, as seen in 
the particle spectra in Figs.[TJ|3] One can encounter the situa- 
tion where some of these states are very weakly populated in 
the reaction, but strongly fed through decay from above-lying 
states. This means that some higher-order 7 rays are not fully 
subtracted in the first-generation method, giving an erroneous 
primary 7-ray spectrum for low Ey. On the other hand, the re- 
action might populate very strongly some low-lying states that 
are more moderately populated by decay from above. One can 
then subtract too much of the 7 rays from these states. 

The latter case is demonstrated for ^"V in the right panel of 



Fig. |4] A state at excitation energy 910 keV with spin/parity 
7+ decaying 100% to the ground state, is strongly populated 
in the neutron pick-up reaction, which favors high-^ transfer 
(here £ = 3, see Ref. [36]). However, it is not so strongly 
populated by 7 decay from above, and the result is that there 
is a vertical valley with zero counts in the primary 7-ray matrix 
at this 7 energy. 

Furthermore, it is important that the populated spin distribu- 
tion is (at least approximately) independent on the excitation 
energy. Else, the bins at high excitation energy will contain 
decay from states with higher spin than the bins at lower ex- 
citation energy, again disturbing the low-energy part of the 
primary 7-ray spectra. For the ^^^Dy(^He, a)'^^Dy reaction, 
the spin population has been extensively studied in Ref. ||23]| 
and references therein. In this specific case, the spin distribu- 
tion turned out to be approximately constant in the excitation- 
energy region investigated. Also, for the lighter nuclei it is 
observed that indeed, the spins are populated with the same 
relative intensities within the error bars. This is seen for ^^V 
in the left panel of Fig.|4j the relative feeding to the low-lying 
states (vertical lines) is approximately the same for the whole 
quasi-continuum region. This is also the case for ^^Ti ifTSll . 

Another potential problem could arise from the finite de- 
tector resolution. To illustrate this, consider a first-generation 
7-ray of 8 MeV, decaying from £, = 10 MeV to Ef = 2 MeV. 
This 7-ray would typically have a resolution of « 250 keV, 
while the particle resolution could be w 150 keV (for 15- 
MeV protons). This means that the weighting function (see 



Sec. II B I for Ef = 2 MeV is about 100 keV broader than 
the excitation-energy resolution at this point. If the situa- 
tion is that there is only one level within Ef, one could then 
"miss" some of the weighting function because it is broader 
than the particle peak. It could also be that the opposite sit- 
uation applies: low-energy 7 rays with, say, 50-keV resolu- 
tion might decay to excitation energies with resolution rang- 
ing from 150 — 300 keV, leading to too narrow weighting func- 
tions. 

We have tested the effects of different resolutions by em- 
ploying a very simple, artificial decay scheme, see Fig. 14 
A hypothetical nucleus with three excited states at Ei — 3.5 
MeV, £2 = 6 MeV, and £3=8 MeV, was assumed to have the 
following decay scheme: 

• from £3: 30% 7„, 20% Yh, and 50% 7.. 

• from £2: 67% 7/, and 33% Ye- 

• from^i: 100% Yf- 

The 7-ray energies involved are: 7, = 2.0 MeV, 7, = 4.5 MeV, 
Yc = 8.0 MeV, Yd = 2.5 MeV, Ye = 6.0 MeV, and Yf = 3.5 
MeV. The first-generation 7 rays from £3 are then 7,, 7^, and 
7c, from E2 Yd and Ye, and from Ei Yf- 

Applying no smoothing for all excitation energies, i.e., 
the 7-ray peaks are 8 functions, the exact result is obtained 
from the first-generation method. Then, we assume 200- 
keV resolution for all excitation energies, but with an energy- 
dependent resolution of the 7-ray spectrum with 50-keV reso- 
lution for 1-MeV 7 rays and 300-keV resolution for 9-MeV 7 
rays, similar to the experimental conditions. This constructed 
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FIG. 14: (Color online). Hypothetical decay scheme of an artificial 
nucleus (see text). 



matrix is shown in Fig. [15] When applying the first-generation 
method on this smoothed matrix, we get the result shown in 



the lower part of Fig. 15 It is seen that it is not exact any more, 
but the differences are small. For example, for £3=8 MeV, 
7j is not a primary transition and should have been completely 
gone in the first-generation spectrum, but still about 5% of the 
counts in the original peak is present. For the and peaks, 
the situation is the same; also so for 7/ at £2 = 6 MeV. This 
means that one might expect leftovers of higher-generation 7 
rays of the order of 5% in the primary spectra. Compared 
to values of experimental errors, which are typically within 
5 — 30%, this is a relatively small effect (the error propagation 
is discussed in detail in Ref. Il25l ). 

To check more thoroughly what effect possible errors in 
the first-generation method might have on the final results, 
namely the extracted level density and strength function, we 
have performed simulations with the generalized version of 
DICEBOX m, as already discussed biefly in Sec. [mA 
Again, we have considered an artificial nucleus resembling 
^^Fe, with a spin distribution of the initial excitation-energy 
levels of 1/2 < y < 13/2, and bin size of 120 keV. Note also 
that equally many negative- and positive-parity states are as- 
sumed above the region of known, discrete levels. 

First, the simulated spectra were folded with the CACTUS 
response functions, and also a Gaussian smoothing was ap- 
plied giving a full width at half maximum (FWHM) of 250 
keV for all excitation energies. These spectra were thus made 
to be as similar as possible to experimental spectra. Then, 
we applied the unfolding technique and the first-generation 
method to obtain the first-generation spectra. We then com- 
pared the extracted first-generation matrix with the true first- 
generation matrix from the simulations, after smoothing the 
true first-generation spectra with a resolution similar to the 
experimental one. Examples of two such matrices for one 
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FIG. 15: (Color online). Matrix of energy levels and 7 transitions in 
an artificial nucleus. Top: original matrix; bottom: result after the 
first-generation method is applied. 



nuclear realization are shown in Fig. 16 The overall good 
similarity between the two matrices is gratifying. However, 
in the low-Ey region, there are significant differences between 
the results extracted from the experiment-like matrix and the 
true, smoothed matrix. In particular, one can see that there 
are some vertical Unes in the experiment-like first-generation 
matrix, e.g., for Ey wl020 keV, that are not present in the 
true matrix. This particular vertical ridge originates from the 
7/2^ state at 1007 keV, which for this nuclear realization is 
strongly populated in the decay cascades at high excitation en- 
ergy. However, at low excitation energies, which corresponds 
to the population from the direct reaction, this state is only 
moderately populated. Thus, its 7 decay is not correctly sub- 
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FIG. 16: (Color online). Simulated first-generation matrix for ^^Fe. 
Top: spectra from folded data set. Bottom: true spectra with smooth- 
ing similar to the experimental resolution. 
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FIG. 17: (Color online). Results from simulated spectra on ^^Fe. 
Top: level density. Bottom: y-ray strength function. 



traded in the first-generation procedure. It is therefore impor- 
tant to exclude such leftovers from higher-generation /rays in 



the further analysis, as mentioned in Sec. 11 C 



We also tested the case with a Gaussian smoothing on the 
particle resolution, and the ideal response on the 7-detection 
side. The extracted level densities and strength functions for 
all three cases are displayed in Fig. [17] In general, the results 
agree very well. We note that there are larger fluctuations in 
the strength function extracted from the experiment-like ma- 
trix as well as the case where only the particle resolution is 
applied, especially in the region below Ey w 2.5 MeV. These 
fluctuations are related to uncertainties in the first-generation 
subtraction procedure and could be due to small variations 



in the shape of the 7 spectra. However, by compressing 
the experiment-like 7 spectra by a factor of two (bin size of 
240 keV), we obtain practically the same shape of the 7-ray 
strength function as from the true first-generation spectra (see 
Fig.p. 

Finally, we have tested the two normalization options (sin- 
gles or multiplicity) of the first-generation method on experi- 
mental spectra to investigate the effect on the extracted level 
density and 7-ray strength function. The result for ^^V is 
shown in Fig. 



18 



We observe that the two options give very 
similar results, only a very few data points are outside the ex- 
perimental error bars. It is clear that the two methods do not 
give any difference in the overall shape of neither the level 
density nor the strength function. 



13 



10" 



> 



LU 

Q. 

ID 
C 

a 

73 

> 



10' 



10 



multiplicity norm, 
a singles norm. 



51i 



2 4 6 8 10 

Excitation energy E (MeV) 



10' 



> 

0) 



realization 1 
realization 2 
Input RSF, no T dep. 




y-ray energy E (MeV) 



o 



u 

1 10 = 



c 



V) 
0] 



:io-' 



multiplicity norm, 
a singles norm. 



i 



51i 



4 6 8 
y-ray energy E (MeV) 



10 



FIG. 18: (Color online). Test of the two normalization options in the 
first-generation method for ^'V. Top: level density. Bottom: y-ray 
strength function. The data are taken from the experiment presented 
inRef |13J. 



C. The Brink hypothesis 

The 7-ray transmission coefficient ^{Ey) in Eq. (fTol is as- 
sumed to be independent of excitation energy (andmus nu- 
clear temperature) according to the generalized Brink hypoth- 
esis f24], as discussed in Sec. |11C| This hypothesis is vio- 
lated when high temperatures and/or spins are involved in the 
nuclear reactions, as shown for giant dipole excitations (see 
Ref. [37] and references therein). However, since both the 
temperature reached and the spins populated are rather low for 
the Oslo experiments, these dependencies are usually assumed 
to be of minor importance in the relatively low excitation- 



FIG. 19: Results from simulated spectra of the artificial nucleus 
'^^Dy. The extracted y-ray strength functions are shown as black 
points and open circles for two different nuclear realizations, while 
the input y-ray strength model is shown as a solid line. 



energy region considered here. 

The effect of the Brink hypothesis has been tested by an- 
alyzing simulated spectra on an artificial nucleus resembling 
'^^Dy. The simulations were again performed with the DICE- 
BOX code 1 34 1 for a specific spin range on the initial excited 
levels, 1/2 <7< 13/2. 

As a first step, a temperature-independent model for the y- 
ray strength function was used as input for the simulations. 
The extracted and input y-ray strength function are shown in 



Fig. 19 As expected, the Oslo method works very well in this 
case. 

In the next test, a temperature-dependent input y-ray 
strength function was used, with temperature T oc y/Ef. In 
principle, it is not possible to disentangle the input level den- 
sity and y-ray strength function anymore, since now we have 



(28) 



P{E,Ey)^p{Ei),^{EuEy). 



This we keep in mind when we use the procedure of Ref. ||25 
in order to extract the level density and y-ray strength func 
tion. 

The extracted y-ray strength functions for two different nu 



clear realizations are shown in Fig. 20 within the excitation- 
energy range 2.1 < E < 6.2 MeV. It is seen that the extracted 
y-ray strength function lies in between the two extremes of 
the temperature-dependent input model, and thus an average 
strength function for the excitation energy region under study 
is found. The shape of the extracted y-ray strength function 
is therefore quite reasonable, although it is clear that the low- 
energy part with 1 <Ey<2.5 MeV must necessarily be quite 
different from the two extremes of the input. 
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FIG. 20: Results from simulated spectra of '^■'Dy. The extracted y- 
ray strength functions are displayed as black and open squares for 
two different nuclear realizations. The input y-ray strength function 
models for Tf = (solid line) and Tj °= ^ (Sn — Ey) (dashed-dotted 
line) are also shown. 



FIG. 21: Results from simulated spectra of the artificial nucleus 
'^^Dy. The extracted y-ray strength functions for two excitation- 
energy regions are displayed as black and open squares. The in- 
put y-ray strength function models for Tj- — (solid line) and Tj °= 
^y {S„ — Ey) (dashed-dotted line) are also shown. 



This is further illustrated in Fig. [2T| where we have ex- 
tracted the y-ray strength function for two separate excitation- 
energy regions: 2.1 < £ < 4.1 MeV and 4.1 < £ < 6.2 MeV. 
It is easily seen that the y-ray strength function for y energies 
between « 1 — 2 MeV is different in the two cases; the higher 
excitation energies lead to a higher temperature of the final 
states and thus a higher strength function. 

However, one should keep in mind that the experimental y- 
ray strength functions have been tested against the assumption 
of temperature dependence for many nuclei, e.g., ^^Sc lfT4l . 
56'57Fe [|10J, '^^■'^^Mojlll and '^^Sn |8|. This is also shown 



asBl 



for '^^Dy 15] in Fig.|22j where the y-ray strength function has 
been extracted for three sets of initial excitation energies. As 
can be seen from the figure, the similarity of the three y-ray 
strength functions is striking. There is, therefore, no experi- 
mental evidence in this excitation-energy region for a strong 
temperature dependence in the strength function. Hence, the 
Brink hypothesis seems to be valid here. 



D. The parity distribution 

As mentioned previously, for both the normalization of the 
level density and the y-ray transmission coefficient, the as- 
sumption of equally many levels with positive and negative 
parity is used. We will in the following investigate this as- 
sumption in detail. 

Using p+ and p_ to denote the level density with positive 
and negative parity levels, the parity asymmetry a is defined 



a = 



P+-P- 
P++P- 



(29) 



which gives — 1 and 1 for only negative and positive parities, 
respectively, and when both parities are equally represented. 

Another expression widely used in literature is the ratio 
p_/p+, which relates to a by 



P+ 



a 
a' 



(30) 



We have considered theoretical parity distributions from 
combinatorial plus Hartree-Fock-Bogoliubov calculations of 
spin- and parity-dependent level densities |39|. Applying 
the definition in Eq. pO] ), we find the calculated parity dis- 
tributions for several Fe, Mo, and Dy isotopes as shown in 



Figs. 23 25 As one might expect from the fact that more or- 
bits are accessible at increasing excitation energy, the parity 
distributions are seen to approach zero as the excitation en- 
ergy increases. However, it is clear that for the lighter nuclei, 
in particular the Fe isotopes, the assumption of zero parity 
asymmetry is not fulfilled in the calculations for excitation 
energies below w 10 MeV. 

We have also looked at other theoretical work such as 
shell-model Monte Carlo calculations |4()J and macroscopic- 
microscopic calculations BTI . In Fig. 2 of Ref. |j41 1, the ratio 
p_ /p+ is shown for ^^Fe, indicating a value of p_ /p+ ~ 0. 1 
at 10 MeV excitation energy. From Fig. 4 in Ref. Ii40i . the 
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FIG. 22: Experimental y-ray strength function of '^^Dy for three sets 
of initial excitation energies. The data are taken from the experiment 
presented in Ref. [5]. 



ratio P-/p+ ~ 0.2 for £ = 10 MeV. In contrast to this, the 
combinatorial calculations of Ref. |39| give p_/p+ ~ 0.7 at 
£ = 10 MeV, which is also in accordance with other mi- 
croscopic calculations based on the Nilsson model and BCS 
quasi-particles ifTTl . where p_/p+ ~ 0.5. These results in- 
dicate considerably more negative-parity states in ^^Fe than 
found in Refs. |40, 41 1. In this specific case, the amount of 
positive-parity states is very sensitive to the position of the 
gij 12 orbital relative to the Fermi level. 

To our understanding, there are currently no experimen- 
tal data on the parity distribution in ^^Fe. However, recent 
measurements on level densities of ~ 2+ and 2^ states 
show no indication of a significantly 



in 5**Ni and ''^Zr 



90^ 



larger amount of states with one of the parities in none of the 
nuclei under study at £ ~ 10 MeV. Also, from the study of 
proton resonances in ^^Sc ll38l . equally many 1 /2+ and 1 /2^ 
states were found, again at £ « 10 MeV. Thus it seems reason- 
able to assume that the parity asymmetry is at least very small 
for these excitation energies, in support of the assumption of 



equal parity as described in Sec. II D 

We would nevertheless like to investigate the impact of the 
assumption of parity symmetry on the calculations of p{S„). 
Let us assume that the spin- and paiity-projected level density 
p{E,J,n) can be described by |41| 



p(£,7,7r)=p(£).g(£,y).^(£,7r), 



(31) 



where p(£) is the total level density at excitation energy £, 
g{E,J) is the spin distribution given by Eq. (27 1, and ^{E,7z) 
is the parity projection factor. According to Eq. ( [T9| , we get 



1 

Do 



-p{Sn)-g{Sn,J^It±l/2)-^{Sn.n,) (32) 



a 
>. 

i: 0.5 

<u 

E 

I, 
<n 

(0 

>. 



-0.5 



(0 
Q. 



a 

^ 0.5 

0) 

E 

I, 

(0 

n 
>. 



-0.5 



(0 
Q. 



>. 

^ 0.5 
E 

I 
n 

.■e -0.5 

(0 

=^ -1 



2 4 6 8 10 12 14 

E (MeV) 

FIG. 23: (Color online). Calculated parity distributions as a function 
of excitation energy calculated for ^^^^'^Fe, from Ref. 1391 . 



for the neutron resonance spacing at S„ reaching states with 
parity ;r, • (— 1)^ = Tit for s-wave neutrons having ^ = 0. Now, 
we define the parity projection factor for positive and negative 
parities as 



and 

using 

Further, 
1 

A) 



^{E,n = n^) 



p+ _l + a 
~P^~' 

p_ _ 1 - a 



= 1. 



(33) 



(34) 



(35) 



p(5„) [g{S„,J = I, + 1/2) +g{S,,J ^I,- 1/2)] ^(5„,;r) 

(36) 
l±a 



: p(5„) [g{S„,J - /, + 1/2) +g{S„,J = /, - 1/2)] 



2 ' 
(37) 
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FIG. 24: (Color online). Same as Fig.|23jfor ''^"'^'^Mo 



FIG. 25: (Color online). Same as Fig.|23]for ^^^-^^Dy. 



which gives 



PiS„) 



(l±a)/2 



Do (/, + l)exp[-(/, + 1)2/2(J2] +/,exp [-I?/2a^] ' 

(38) 

using Eq. ( [27] l. If the target nucleus in the neutron-capture 
experiment has positive parity in the ground state, the factor 
( 1 + a)/2 is used, and for negative ground-state parity we use 
(l-a)/2. 



For several key cases, we have applied Eq. (38 i for calcu- 
lating p{S„) and compared to the result using Eq. (20i. For 
example, for ^^Fe with neutron resonance spacing Do = 6.5 
keV at S„ = 10.044 MeV IIS, and using the spin cutoff pa- 
rameter g{S„) ~ 3.93 from the prescription of Ref. |29|, we 
obtain p(5'„) =2518 MeV^' if the assumption of equal parity 
is used, and p{S„) = 2939 MeV^' if we correct for the parity 
asymmetry a = 0.14 predicted by the combinatorial model of 
|!39|. Thus, including the parity asymmetry gives about 17% 
higher level density at S„. For '^''Mo, with Dn = 105 eV ll33l. 
a{S„) = 5.15, and a = -0.017 |39|, we get p(5„) = 1.01 • 10^ 
MeV^^ when no parity asymmetry is taken into account, and 
p{S„) = 1.03 • 10^ MeV^' when the parity asymmetry is con- 
sidered; only a change of sa 2%. 



To check the effect of the parity asymmetry on the nor- 
malization, we have renormalized the data on ^^Mo using 
the above-mentioned values for p{S„) with and without par- 
ity correction. Since the estimated a was so small for this 
case, we have also assumed a much larger parity asymmetry 
of a = -0.5 leading to p(5„) = 2.02 • 10^ MeV-', a factor of 
2 larger level density than the one assuming equal parity. As- 
suming that a ~ 0.5 gives p(5„) = 6.74- 10^ MeV"', roughly 
a factor of 2/3 reduction compared to the parity-symmetry 
case. 

The resulting level density and strength function are shown 
in Figs. 26 Here, it is easily seen that a small parity asym- 



metry does not give any significant changes of the normal- 
ization in neither the level density nor the strength function. 
However, for a — —0.5 corresponding to a factor of 2 larger 
p{Sn) gives an overall larger level density for excitation en- 
ergies larger than ss 3 MeV, and the steeper slope is reflected 
also in the strength function. In addition, one sees a suppres- 
sion in the y-ray strength function for Ey <6 MeV, a direct 
consequence of the changes in the level density. The oppo- 
site is true for a — 0.5; here, p{S„) and thus the slope of the 
level density is reduced. Consequently, the slope is reduced 
also in the y-ray strength function, but the absolute value for 
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FIG. 26: (Color online). Effect of parity asymmetry on the ex- 
perimental level density (upper panel) and strength function (lower 
panel) of ^^Mo. The data are taken from the experiment presented in 
Ref. 112J. 



Ey<6 MeV is in this case increased with respect to the parity- 
symmetry case. Note that such large values of a represent ex- 
treme cases; both from experimental data |'38','425 and several 
theoretical estimates (e.g, |11, 39|) the asymmetry around 5„ 
should be smaller than typically ±0.2. 

The assumption of parity symmetry influences also the nor- 
malization of the 7-ray transmission coefficient. To take into 
account the parity distribution, one can modify Eq. ([25]l ac- 



p{E-Ey,JuTti)=p{E-Ey)-^{E-Ey,Ki)-g{E-Ey,Ji). 

(39) 

We will restrict ourselves to consider t= \ radiation only {E 1 
and Ml), since in general, this multipolarity is expected to 
give by far the largest contribution to the strength function in 
the quasicontinuum region (see, e.g., studies by Kopecky and 
Uhl 1.321). 

For E\ radiation, the parity of the final state is opposite to 
the initial state. In s-wave neutron capture experiments, the 
parity of the target nucleus' ground state is equal to the parity 
of the neutron resonances of the created compound nucleus. 
Therefore, the accessible parity of the final states must be the 
opposite of the initial state. For Ml radiation, the accessible 
final states must have the same parity as the initial state. 

Based on Eq. ( [2T| i and taking the parity distribution into 
account, one finds 



{Ty{S„J,±\/2,n,)) 



B 



£y=0 



27tp{S,„Ir±l/2,n,) 

dEy [^El (Ey) ^_/+ + =Tm1 {Ey) ^+/_] 



xp{S„-Ey) ^ g{S„-Ey,Ir±l/2+J). (40) 

j=-i 

for target nuclei with positive/negative ground-state parity. It 
is easily seen that if ^+ — — Eq. (26 1 is restored. 

The Oslo method does not enable the separation of ^^{Ey) 
into its £1 and Ml components. We have therefore applied 
models for the level density and the £1 and Ml strength func- 
tions in order to investigate the influence of including par- 
ity on the integral in E q. ([40| . For the level density, we 
have used the results of M39I . For the £1 strength function, 
we have used the the Kadmenskii, Markushev and Furman 
(KMF) model |43 1 with a constant temperature Tf = 0.3 MeV, 
while for the Ml component (the spin-flip resonance ll22l ') we 
have applied a Lorentzian shape (see Ref. |33|). We use again 
^^Mo as a test case, with experimental GDR parameters taken 
from Ref. Il33l and Ml Lorentzian parameters from systemat- 
ics ||33]| . The target nucleus ^^Mo in the («,7)^''Mo reaction 
has spin/parity 5/2+. 

First, we take the parity distribution from Ref. Il39l as 
shown in Fig. |24]and apply in the integral of Eq. ( [40| . The 
resulting value is only about 1% smaller than the value ob- 
tained assuming — = 1/2 for all excitation energies, 
implying that the effect of parity is negligible in this case. As 
a further test we used the parity distribution of ^''Fe shown 
in Fig. 23 on ^^Mo; then, we found a w 47% reduction using 
Eq. (40 1 compared to the parity-symmetry case. For the very 
extreme cases assuming that a = 1 for all E (allowing only 
for positive-parity states and Ml transitions), we get about a 
factor of 3 reduction, and using a = — \ for all E (allowing 
only for negative-parity states and E 1 transitions), we obtain 



about 67% increase of the normalization integral in Eq. (40 1 
relative to the case where no parity distribution is considered. 

To summarize, for nuclei in the Mo mass region and above, 
we find only small corrections of the order of a few percent 



18 



x10"- 



40 



^ 35 



> 

c 

u 
c 



30 



O) 

c 

0) 



extracted f.g. spectra, equal n 
true f.g. spectra, equal % 
extracted f.g. spectra, asym. % 
true f.g. spectra, asym. n 
input RSF 




y-ray energy (MeV) 

FIG. 27: (Color online). Extracted strength functions from simulated 
data for an artificial Fe nucleus, with and without parity asymme- 
try in the level density. 



to the estimation of p(5„) and also for the absolute normal- 
ization of the strength function when using realistic parity 
distributions (compared to the extreme cases). However, for 
light nuclei such as Fe, effects of the parity distribution could 
be significant on both the level-density and strength-function 
normalization, of the order of 30-50%. 

The parity distribution could also potentially influence the 
y-ray strength function in other ways than just the normal- 
ization. In the quasi-continuum region, one expects that £1 
transitions largely dominate as long as the parity asymmetry 
is close to zero. However, if the parity asymmetry is large. 
Ml transitions will be favored over those of £1 type. We 
have investigated this using DICEBOX |34| to generate simu- 
lated data on an aitificial nucleus resembling ^^Fe, as already 
described in Sec. 



IllB 



We have applied a symmetric parity 
distribution in the first case, and an asymmetric parity distri- 
bution in the second case. Note that the parity distribution is 
implemented directly in the level density only; the input 7-ray 
strength function is kept fixed. 

The resulting 7-ray strength functions are shown in Fig. 27 
Visually, it is hard to see any difference between the 7-ray 
strength functions with and without parity asymmetry. By tak- 
ing the average 7-ray strength function for 7 energies between 
1 .4 — 3 MeV, we find an average increase of about 20% for the 
case with parity asymmetry relative to the equal parity case, 
using the true first-generation spectra. The statistical fluctu- 
ations here are typically less than 10% for Ey = 2 — 3 MeV, 
and up to 40% for 7 rays below 2 MeV. If we use the extracted 
first-generation spectra from the folded data, the difference 
between the two cases is approximately 6-7%, less than the 
statistical fluctuations which are larger than when using the 
true first-generation spectra. For higher Ey, there is no effect 



of parity within the fluctuations. We therefore conclude that it 
could be an effect of w 20% on the low-energy part of the 7- 
ray strength function. Hence, we find it reasonable to believe 
that even a considerable parity asymmetry on the level density 
as in the ^^Fe case will not drastically change the shape of the 
7-ray strength function. 



In Fig. 27 we note that there is an enhanced strength for 
7-ray energies below w 3.5 MeV compared to the input 7-ray 
strength function, regardless of the parity distribution. This 
is in fact due to the spin range of the initial levels, which is 
further investigated in the following section. 



E. The spin distribution 

The uncertainty in the spin distribution may lead to er- 
rors in the Oslo method in three ways: (;) the extraction of 
first-generation 7 rays and subsequent effects on the extracted 
level density and 7-ray strength function, (//) the estimation of 
p{S„), and (///) the spin range accessed experimentally (typi- 
cally — 8/j depending on the reaction and target spin) com- 
pared to the true, total spin distribution. These issues will be 
discussed in the following. 

We have once more relied on simulated data applying the 
DICEBOX code I34il in order to test the sensitivity on the 
final results with respect to the spin range of the initial lev- 
els. We simulated again a light nucleus that resembles ^^Fe, 
as we expect that light nuclei would be most sensitive to the 
initial spin population. This is because the higher spin levels 
are missing at low excitations. Again, the critical energy was 
set to 2.2 MeV, as there are no available data in the literature 
about levels with higher spins above this energy. As before, 
all levels had the same probability of population. The input 
7-ray strength function model was independent of excitation 
energy. 

Three spin ranges were used: 1/2 < J < 1 /2, 1/2 < J < 
13/2, and 7 /2 < 7 < 13/2. This means that the highest possi- 
ble spin reached at the final excitation energy was 9/2 or 15/2 
depending on the maximum allowed initial spin. The weight- 
ing of the allowed spins followed the spin distribution of the 
input level density, which in this case was from the micro- 
scopic calculations of Ref. |39|. For all initial spin ranges we 
applied the full Oslo method, and in parallel we extracted the 
level density and 7-ray strength function from the true first- 
generation spectra. The resulting level densities and 7-ray 



strength functions are displayed in Figs. 28 and 29 respec- 
tively. 

As can be seen from Fig.|29] there is an enhanced strength 
at low 7 energies for the spin windows including the highest 
spins. Also, by looking at the average multiplicity (M) at 7.0 
MeV of excitation energy, we get 1.8 for 1/2 < 7 < 7/2, 2.2 
for 1/2 < 7 < 13/2, and 2.5 for 7/2 < 7 < 13/2. At By = 
2.1 MeV, we find that the 7-ray strength function is increased 
with a factor of « 1.7 for the spin range 1/2 < 7 < 13/2, 
and a factor of w 2.6 for the spin range 7/2 < J < 13/2 with 
respect to the input 7-ray strength function. This is not an 
artifact from the unfolding or first-generation method, as it is 
also seen in the extracted 7-ray strength functions from the 
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FIG. 28: (Color online). Extracted level densities from simulated 
data for initial spin range 1/2 < 7 < 7/2 (top panel), 1/2 < 7 < 13/2 
(middle panel), and 7/2 < 7 < 13/2 (bottom panel) on the initial 
levels. 



FIG. 29: (Color online). Extracted strength functions for spin range 
l/2<y<7/2 (top), l/2<y< 13/2 (middle), and7/2<y< 13/2 
(bottom) on the initial levels. 



trae first-generation spectra. 

The explanation for the observed behavior is probably con- 
nected to three issues: (/) the dominance of dipole radiation, 
which carries L — I; (//) the applied spin restriction on the 
initial levels; (///) the low level density in this mass region at 
low excitation energies, especially for high spins. Let us take, 
as an example, the population of a 13/2 level at high excita- 
tion energy. If it is to decay to a low-lying level, involving 
a high-energy 7 transition, a level with appropriate spin must 
be present at this excitation energy. This is not necessarily the 
case as there are only about 2—10 levels per 120 keV for exci- 
tation energies below 4.7 MeV. Among these few levels, there 
are probably no high-spin levels at all. This leads to a higher 
average multiplicity and an enhanced probability for decaying 
with low-energy 7 rays. 



From Fig. 29 it is clear that the enhancement in the ex- 



tracted 7-ray strength function is not present in the input 7-ray 
strength function, but is an effect of the three factors men- 
tioned above. This means that the simple factorization in 



the input 7-ray strength function. The extracted 7-ray strength 
function, which could be considered as an "effective" 7-ray 
strength function, is however fully capable of reproducing the 



true primary 7-ray spectra, and it is seen from Fig. 28 that the 
extracted level density is very reasonable indeed. 

However, it is not obvious whether the spin range of the 
inital states is the full explanation of the observed low-energy 
enhancement in light and medium-mass nuclei. Further inves- 
tigations are therefore needed to clarify this issue. It should 
also be noted that similar tests have been performed for an ar- 
tificial nucleus resembling '^^Dy. No such enhancement was 
seen in this case, in agreement with experimental findings in 
this mass region. This is not surprising since in these nuclei, 
the level density is much higher and relatively high spins are 
available already at low excitation energies. 

The quantity p{S„) is calculated assuming a bell-Iike spin 
distribution according to Ref. Il27l given by Eq. ( 27 1 and us- 



Eq. ( 10 1 of the first-generation spectra is not able to reproduce 



ing a model for the spin cutoff parameter a usually taken from 
Ref. II27II or from Ref |29|. Both these assumptions could in 
principle be a source of uncertainty, as it is hard or even im- 
possible to measure experimentally the total spin distribution 
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at high excitation energy. Thus, the theoretical spin distribu- 
tions are seldom constrained to data for excitation energies 
higher than typically 2 — 3 MeV. 



In Fig. 30 various spin distributions for ^'^Sc are shown, cal- 
culated at an excitation energy of 8.0 MeV. In the two up- 



per panels, the spin distribution given in Eq. ( 27 1 has been 



used, but with the expression for the spin cutoff parameter 
of Refs. lIZTl l28l in panel a) and the formalism of Ref. |,29J 
in panel b). In panel c) the spin distribution of the spin- 
dependent level densities of Ref. p44l are shown. Here, the 
authors also have assumed a bell-shaped spin distribution ac- 
cording to Eqs. (7) and (8) in Ref. Il44l . It is clear from the 
figure that the spin distributions in panel b) and c) are broader 
with centroids shifted to higher spins compared to the one in 
panel a). 

In panel d), the spin distribution of the calculated spin- and 
parity-dependent level density of Ref |45| is shown. There 
are no underlying assumptions for the spin distribution in 
these calculations. It is seen from this distribution that there is 
a significant difference in the relative numbers of states with 
spin and 1 . The normalization procedure for the level den- 
sity described above is especially sensitive to such variations 
at low spin if the neutron resonance spacing Dq is measured 
from a neutron capture reaction where the target nucleus is 
even-even, that is, with zero ground-state spin. Then the states 
reached in neutron capture can only have spin 1 /2+, and the 
number of all other states must be estimated using a certain 
spin cutoff parameter, introducing a larger uncertainty in the 
calculated p{S„). Therefore it is preferred to calculate p{S„) 
from both Do and Di resonance spacings if possible, since in 
the latter, also states with 3 /2^ are reached for target nuclei 
with //^ = 0+, and will therefore decrease this uncertainty. 

We have investigated the relative difference in the value 
of p{Sn) for several nuclei using the spin-cutoff parameter 
of [27] with global parameterization of |28| (referred to as 
Gilbert & Cameron), and the prescription of |29| (hereafter 
called von Egidy & Bucurescu), see Tab.|I] For the cases inves- 
tigated here, the general trend is that the Gilbert & Cameron 
approach leads to a lower value for p{S„) than the von Egidy 
& Bucurescu parameterization. However, this depends on the 
spins reached in the («, 7) reaction. For example, for V, we 
note that the obtained p {S„ ) using the spin cutoff of von Egidy 
& Bucurescu is lower than that of Gilbert & Cameron. This is 
easily understood from the fact that in this case, relatively high 



spins are populated in the (n, 7) reaction. We see from Fig. 31 



that the Gilbert & Cameron parameterization gives lower rel- 
ative values for states with spin 1 1/2 and 13/2 compared to the 
von Egidy & Bucurescu distribution. Thus, one divides by a 
smaller number in Eq. (20 1 and p{S„) gets larger. 



One can conclude that the two approaches studied here give 
a w 10 — 50% change in the resulting p (5„), while other calcu- 
lations f39l indicate even larger variations. Using this HFB- 
plus-combinatorial approach gives typically a factor of two 
change compared to paciSn) for most cases and up to a factor 
of 3.7 for the extreme case of '^^Dy (see Tab. II in Ref. |39J). 
The effect of these differences is of course an uncertain deter- 
mination of p{S„) and thus of the slope of the level density 
and the 7-ray strength function, in a similar fashion as for the 
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FIG. 30: Relative spin distributions calculated for E = 8.0 MeV of 
''^Sc (see text). 



parity dependence (see Sec. HID and Fig. 26 1 



There is also a question of how the spin range accessed ex- 
perimentally relative to the true, total spin range might influ- 
ence structures in the level density. In the analysis, we assume 
that the excluded spins contribute on average with a scaling 
factor of the level density, which is corrected for by normal- 
izing to p{S„), i.e. the total level density for all spins at this 
energy. However, this rehes on the hypothesis that structures 
in the level density due to, e.g., nucleon pair breaking are not 
severely affected by the spin window. 

To test this, we have performed simplistic calculations with 
the code COMBI |46|, which is based on a microscopic, com- 
binatorial model using Nilsson single-particle energies and the 
concept of quasiparticles from the nuclear Bardeen-Cooper- 
Schrieffer (BCS) theory. The chosen spin windows were 
< y < 6 and < 7 < 30 (in units of h). Details of the model 
can be found in Ref. |46l. 

The resulting calculated level densities for ^^Fe and ^^^Dy 



are displayed in Figs. |32j and l33j respectively. For ^''Fe, is 



it clear that most of the states are found within the spin range 
of — 6h. In fact, the relative difference between the large 
and the small spin window is at most « 30%, and it is clear 
that there are no significant structural differences for the two 
cases. One would therefore not expect any severe problems 
with normalizing the level density of such light nuclei to the 
total level density p{S„). However, for the '^^Dy case, the 
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TABLE I: Variation of the calculated p{S„) using different spin cutoff parameters. The target spin in the {n, y) reaction is denoted by If, and 
El is the total backshift for the back-shifted Fermi-gas model while a is the level-density parameter. All level spacings (Do) are taken from 
L33J . The parameter rj is the ratio PEsiSn)/ Pcc{Sn)- 
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FIG. 31: Relative spin distributions calculated at 5„ = 11.051 MeV 
of ^'V (see text). 

deviation of the level densities calculated within the two spin 
ranges can be considerable and is increasing with excitation 
energy, up to a factor of « 2.5 for £ w 7.5 MeV. It is nev- 
ertheless very gratifying that, indeed, the gross structures of 
the level densities are very much alike, so that extracted infor- 
mation from the experimental level density, such as the onset 
of the pair-breaking process, is probably reliable (the larger 
spin window leads to a smoothing of the steps). However, ex- 
perimental work to investigate the spin distribution further is 
highly desirable. 

IV. SUMMARY AND CONCLUSIONS 

In this work we have addressed uncertainties and possi- 
ble systematic errors that one can encounter using the Oslo 
method. The main steps of the method have been outlined, 
and the assumptions behind each step are investigated in de- 



FIG. 32: Calculated level density of Fe for two spin ranges: 
< y < 6 (dashed line) and < 7 < 30 (solid line) compared to 
experimental data (open circles, from Ref. II II ) normalized to the 
total level density at S„. 



tail. Our findings indicate that although the Oslo method is in 
general very robust, the various assumptions it relies on must 
be carefully considered when it is applied in different mass 
regions. The results can be summarized as follows. 

The unfolding procedure. We have investigated the effect of 
changing the total absorption efficiencies of the Nal(Tl) crys- 
tals up to w 20%. We found that the extracted level density is 
hardly affected at all, while the 7-ray strength function gets a 
somewhat different slope. The quantitative effect on the 7-ray 
strength function coiTesponds directly to the imposed changes 
in the efficiencies. 

The first- generation method. In general, the assumption of 
thermalization can be questioned at excitation energies above 
« 10 MeV. Also, variation in the populated spin range as a 
function of excitation energy may lead to an erroneous sub- 
traction, especially of yrast states and other strong transitions 
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FIG. 33: Same as Fig. 32 for '^^Dy, with experimental data from 
Ref. 0. 



at low 7-ray energies. This is particularly pronounced at low 
excitation and low 7-ray energies, where the direct reaction 
may favor or suppress the population of some levels compared 
to feeding from 7 decay. In addition, the finite, and for some 
energies, the mismatch of the detector resolution of the parti- 
cle telescopes and CACTUS might in principle lead to prob- 
lems with determining correct weighting functions. However, 
tests on simulated spectra show that this effect is probably of 
minor importance. The two different normalization options in 
the method turn out to give very similar results. 

The Brink hypothesis. Since the extraction of level density 
and 7-ray strength relies on the Brink hypothesis, one would 
expect problems if there was an effective temperature depen- 
dence in the 7-ray strength function. From the results of sim- 
ulated data, it is clear that especially the region of low Ey is 
affected when the 7-ray strength function is temperature de- 
pendent, and the extracted strength is found to found to lie in 
between the two temperature extremes considered. However, 
we note that there is, so far, no experimental evidence for any 
strong temperature dependence of the 7-ray strength function 
in the excitation-energy region in question (below ss 10 MeV). 

The parity distribution. The assumption of equally many 
positive- and negative-parity states may break down, espe- 
cially for light nuclei, and would then affect the normaliza- 
tion of both the level density and the 7-ray strength function. 
If the parity asymmetry is large, the value of the normaliza- 
tion point p{S„) might change with a factor of « 2, which 
will consequently change the slope of the level density and 7- 
ray strength function. For heavier nuclei, however, the parity 
distribution at S„ is expected to be quite even, and experimen- 
tal data on lighter nuclei suggest that the parity asymmetry 
should be small also here. Therefore, we find it reasonable to 



believe that the error in p{S„) and on the absolute normaliza- 
tion of the 7-ray strength function due to this assumption does 
not exceed 50%. 

Another issue is the influence on the relative contribution 
of £1 and Ml transitions, as a large parity asymmetry will 
favor Ml transitions. Tests with simulated data indicate, how- 
ever, that the effect is not large even with a considerable parity 
asymmetry. 

The spin distribution. From simulations on light nuclei with 
different spin ranges on the initial levels, it is seen that the 
ranges including higher spins might lead to an enhanced 7- 
ray strength function at low 7-ray energies. This is under- 
stood from considering the relatively low level density and 
the matching of spins between the initial and final levels for 
dipole radiation. This could explain, at least partially, the ob- 
served enhanced low-energy strength in experimental data for 
light and medium-mass nuclei. 

The spin distribution is one of the largest uncertainties in 
the normalization, since the determination of the slope of the 
level density and 7-ray strength function strongly depends on 
the relative intensities of the populated spins at 5„. Assuming 
a bell-shaped spin distribution with various global parameter- 
izations of the spin cutoff parameter give up to « 50% change 
in p{S„); however, microscopic calculations indicate that this 
value can vary within a factor of 2 or more, especially for 
nuclei in the rare-earth region. This will also consequently 
influence the 7-ray strength function. 

The experimentally reached spin window in OCL experi- 
ments is typically — 8/i, which means that the extracted level- 
density data in fact represent this narrow spin range. This is 
usually compensated for by normalizing to the total level den- 
sity at Sn ■ This is not a big effect for light nuclei with relatively 
few high-spin states. Even in the case of rare-earth nuclei, our 
calculations indicate that this scaling works well as the main 
structures in the level density are indeed present in the levels 
of a rather narrow spin range. 
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